Contact values of the radial distribution functions of additive hard-sphere mixtures in 

d dimensions: A new proposal 
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The contact values gij{(Tij) of the radial distribution functions of a d- dimensional mixture of 
(additive) hard spheres are considered. A 'universality' assumption is put forward, according to 
which Qiiioij) = G{ri,Zij), where G is a common function for all the mixtures of the same di- 
mensionality, regardless of the number of components, rj is the packing fraction of the mixture, 
and Zij = {aiCTj/aij)(a'^~^)/{(j'^) is a dimensionless parameter, (ct") being the n-th moment of the 
diameter distribution. For d = 3, this universality assumption holds for the contact values of the 
Percus-Yevick approximation, the Scaled Particle Theory, and, consequently, the Boublfk-Grundke- 
Henderson-Lee-Levesque approximation. Known exact consistency conditions are used to express 
G(r], 0), G{r], 1), and Gijj, 2) in terms of the radial distribution at contact of the one-component sys- 
tem. Two specific proposals consistent with the above conditions (a quadratic form and a rational 
form) are made for the ^-dependence of G{ri, z). For one- dimensional systems, the proposals for the 
contact values reduce to the exact result. Good agreement between the predictions of the proposals 
and available numerical results is found for ci = 2, 3, 4, and 5. 



I. INTRODUCTION 

It is well known that there exists a close connection 
between the thermodynamic and structural properties of 
classical fluids. In fact, for pair-wise additive intermolec- 
ular potentials, all the thermodynamic functions may be 
expressed in terms of the radial distribution functions 
(rdf). The expressions are particularly simple for hard- 
core fluids, since in that case the internal energy reduces 
to that of the ideal gas and in the pressure equation it 
is only the contact values rather than the full rdf which 
appear explicitly. Therefore, knowledge of the contact 
values of the rdf in hard-core fluids, which we will denote 
by gij{(Tij) (where in general i and j label species and 
tJij is the distance of separation at contact between the 
centers of two interacting fluid particles, one of species i 
and the other of species j), suffices to obtain the equation 
of state (EOS) of these systems. In the case of a single 
component hard-core fluid, if the EOS were known then 
it would be straightforward to infer the contact value of 
the radial distribution function. In contrast, if one were 
given the EOS of a multicomponent hard-core mixture 
such contact values could not be determined in a unique 
way. Up to this day, no exact expressions neither for 
the contact values of the rdf nor for the EOS (except for 
the case of one-dimensional systems, i.e., hard rods) are 
known, although various approximate theories, empirical 
efforts, and computer simulations have been carried out 
in connection with this problem. One-component sys- 
tems are of course easier to handle and this has meant 
that studies for mixtures are much scarcer. Perhaps the 
most successful theoretical (analytical) approach to this 
issue in the case of additive hard-sphere mixtures (d = 3) 
is the exact solution of the Percua-rYevick (PY) equation 
carried out by Lebowitz in 1964.lil This analytical solu- 



3en recently proposed by Hemder- 
'^^ Matyushov and Ladanyijl3 and 



tion, which among other things yields explicit expressions 
for the contact values of the rdf, as well as for the virial 
and the compressibility routes to the EOS, is at the ba- 
sis of the celebrated (and empirically derived) Boublikp; 
Mansoori-Carnahan-Starhng-Leland (BMCSL) EOS,!-! 
considered to be a rather accurate EOS for hard-sphere 
mixtures. In his paper ,□ Boublfk also introduced an ap- 
proximation for the contact values of the rdf (in fact an 
interpolation between the PYj-ijcsults and the ones of the 
Scaled Particle Theory (SPT)ala) that later was indepen- 
dently proposed by Grundke and HendersonQ and Lee 
and Levesque.Q This approximation, which we will re- 
fer to as BGHLL, leads precisely to the BMCSL EOS 
when substituted into the statistical mechanical formula 
for the pressure equation. Refinements of the BGHLL 
approximationrjlja 
son and Chan,B'E[L— . 

Barrio and SolanaM to cope with some deficiencies of the 
BMCSL EOS in the so-called colloidal limit of binary 
hard-sphere mixtures. 

As far as we are aware, there are no reported approx- 
imations for_ftj(cry) with c? 7^ 3, except that of Jenkins 
and [Mancinill3 in the case of hard-disk mixtures and our 
ownta for c?-dimensional mixtures. The latter approxima- 
tion, however, was introduced only as a means to derive 
a proposal for the EOS of mixtures. In fact, while this 
EOS pres ents .an exc elleut agreeme^it with simulations 
for d = 2jmi d = sjmm d = 4,E2I and^^ 5,E2l the 
expressions for gij{crij) are less accurate£3Ej It is the 
major aim of this paper to propose new (improved) ap- 
proximate expressions for gij(aij), for arbitrary mixtures 
and arbitrary dimensionality, that, apart from satisfy- 
ing known consistency conditions, are sufficiently gen- 
eral and flexible to accommodate any given EOS for the 
single fluid. A key aspect of the present approach, also 
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included in our previous work,£3 is that we will assume 
from the very beginning a kind of universal behavior of 
the contact values whicfuat least holds also for the so" 
tion of the PY equation,EI for the[S|PfT approximation, 
and for the BGHLL interpolationaalD in the case of mix- 
tures of hard spheres {d ~ 3). This means that, once 
the dimensionality and the packing fraction are fixed, 
the expression for the contact values gij {uij ) for all pairs 
ij is the same, irrespective of the composition and the 
number of components in the mixture. This expression 
must comply with three consistency conditions related 
to the point particle, equal size, and the colloidal limits, 
respectively. Two functional forms (a quadratic one and 
a rational one) which are sufficiently representative will 
be examined. Their merits will be assessed from a com- 
parison with available simulation data as well as with 
respect to the performance of the EOS obtained from 
them. In this latter issue, we will show that a paradoxi- 
cal result is obtained. What we find is that, contrary to 
what one could possibly expect, better contact values do 
not necessarily mean more accurate EOS and that even 
two different expressions for gij (aij ) may lead to exactly 
the same EOS. 

The paper is organized as follows. In Sec. II we re- 
call the known consistency conditions and introduce the 
new proposals for the contact values of the rdf. Section 



III deals with the comparison between our contact val- 
ues and ensuing EOS and simulation results. We close 
the paper in Sec. IV with further discussion and some 
concluding remarks. 



II. THE PROPOSAL 

Let us consider a mixture of hard spheres in d dimen- 
sions with an arbitrary number of components. The hard 
core of the interaction between a sphere of species i and a 
sphere of species j is cry = ^{<Ti + aj), where the diameter 
of a sphere of species i is an — ai . The number density 
of the mixture is p and the mole fraction of species i is 
Xi = Pi/ P- From these quantities one can define the pack- 
ing fraction rj = Vdpicr'^), where Vd = (n/Ay/^ /r{l + d/2) 
is the volume of a d-dimensional sphere of unit diameter 
and (cr") = 2:^(7" denotes the moments of the diame- 
ter distribution. 

In a hard-sphere mixture, the knowledge of the contact 
values gij{<Jij) is important for a number of reasons. For 
example and as stated above, the availability of gij{(Jij) 
is sufficient to get the equation of state (EOS) of the 
mixture via the virial expression 



the entire rdf gij{rXJa\ the Generalized Mean Spher- 
ical ApproximationEll and in the Rational Function 
Approximation.c3 In a different context, they are impor- 
tant as well to implement the Enskog . kin o tic theory both 
for elastic and inelastic hard spheres. t3E3 

The exact form of 5y (cry ) as functions of the packing 
fraction, the set of diameters {o'k}^ and the set of mole 
fractions {xk\ is only known in the one-dimensional case, 
where one simply has gij{(Jij) = (1 — ?y)~^. Consequently, 
for d > 2 one has to resort to approximate theories or em- 
pirical expressions. From that point of view, it is useful 
to make use of exact limit results that can help one in 
the construction of approximate expressions for gijicij)- 
Let us consider first the limit in which one of the species, 
say i, is made of point particles that do not occupy vol- 
ume, i.e., (Ti 0. In that case, gii{(Ji) takes the ideal gas 
value, except that the available volume fraction is 1 — 77: 



lim Qijiaj) = . 

ai^o^^ ' 1-77 



(2) 



An even simpler situation occurs when all the species 
have the same size, {cr/c} — > cr, so that the system be- 
comes equivalent to a one-component system. Therefore, 



lim 5y(cry) 



(3) 



i,3 



where (?(cr) is the contact value of the radial distribution 
function in the one-component case. Equations (|^) and 
represent the simplest and most basic conditions that 
gij{(7ij) must satisfy. There[-|iS|-a , rumnb f jr . pf rOther less 
trivial consistency conditions.fflliil'EitifEaEjEa'EII Here we 
consider the condition stemming from a binary mixture 
in which one of the species, say i = 1, is much larger than 
the other one, i.e., (Ti/a2 00, but occupies a negligible 
volume, i.e., xi{ai/(T2Y 0- In that sphere of 

specj^^^is seen as a wall by particles of species 2, so 



lim [512(^12) - 2''-i77.g22((T2)] = 1. (4) 

a:i((Ti/cr2)''^0 

Also in that Hmit,iii'0 Ingu^ai) ^ (7i/a2, but we will 
not make use of this condition here. 

Our purpose now is to propose approximate expres- 
sions for gij [uij ) of hard-core mixtures with an arbitrary 
number of components and arbitrary dimensionality d, 
that satisfy the consistency conditions (^)"(^. First, we 
assume that the dependence of gij{crij) on the parame- 
ters {fTfe} and {xfc} takes place only through the scaled 
quantity Zij = {aiaj /aij){(T'^^^) / (a^) . More specifically. 



where = p/pksT is the compressibility factor 

of the mixture, p being the pressure, the Boltz- 
mann constant, and T the absolute temperature. The 
contact values gijifJij) are also needed to generate 



ff»i(cry) = G{T],z,j), 



(5) 



where the function G{r], z) is universal in the sense that 
it is a common function for all the pairs ij, regardless of 
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the composition and number of components of the mix- 
ture. Of course, the function G{ri, z) is different for each 
dimensionahty d. 

The ratio ^ = {a'^~^) / {a'^) can be understood as a 
'typical' inverse diameter or curvature of the particles of 
the mixture. The parameter z~'^ = {(J^^ + aJ^)/2^ rep- 
resents then the average curvature, in units of ^, of a 
particle of species i and a particle of species j. Accord- 
ing to Eq. (|^), if two different pairs ij in two different 
mixtures (with the same packing fraction) have the same 
dimensionless average curvature z~i^, then they have the 
same contact value of the rdf. i— , 

The ansatz (^) includes the one used by usjtB where 
G(?7, z) was approximated by a linear function of z. A 
particular case of this linear form is the proposaJ made 
by Jenkins and Mancini for hard-disk mixtures :t3 

In the three-dimensional case, Eq. (|^) is alsO|-|Com- 
patible with the soiution of the PY equationjll the 
SPT appjuccdmationpla and, consequently, the BGHLL 
proposaLaEfLl More specifically, G(ry, z) is a linear func- 
tion of z in the PY approximation and a quadratic func- 
tion in the SPT and BGHLL approximations: 



Giv,z) 



1 



1 - r/ 2 (1 - ry)2 



(7) 

where Apy = 0, Aspt = f, and Abghll = |- These 
three approximations are consistent with (H) and (2), 
but only the SPT is also consistent with condition (4). 
The approxiraation referred to as the SPT-virial route 
by Rosenfeldfl adopts also the scaling form (^), namely 
G{ri,z) = (1 — 7y)~^ exp[3z?7/2(l — ry)], but it does not 
comply with condition (^). 

Once we adopt the ansatz (|^), the limits in (||)-(|^) 
provide very useful constraints on the z-dependence of 
G. First, Zii — > in the limit ct; — > 0, so that insertion of 
Eq. (|) into d) yields 

G{v,0) ^ 



1-77 

Next, if all the diameters are equal 



Eq. implies that 



G(77,l)-.g(a). 



(8) 
1, so that 

(9) 



Finally, in the limit considered in Eq. , we have Z22 ^ 
1, 2:12 — > 2. Consequently, 

G(r/,2) = l + 2'^-ir/G(??,l). (10) 

Thus Eqs. (P)-(pO|) provide complete information on the 
function G at 2: = 0, 2 = 1, and 2 = 2, in terms of the 
contact value g{a) of the one-component rdf. 

The proposal made in Ref. |l^ consists of assuming a 
linear dependence of G on 2 that satisfies the require- 
ments d) and (|): 



G{v,z) = 



9i<y)- 



(11) 



If irL-the two-dimensional case we take Henderson's 
valucH g{a) = (1 - 7?7/16)/(l - 77)^, Eq^(fTl) reduces 



to Jenkins and Mancini's approximation ,113 Eq. m). In 



general, Eq. does not satisfy Eq. (|10|). However, the 
ansatz ( |ll| ) was used in Ref. ^ only as a means to ob- 
tain the EOS from Eq. @). The resulting EOS exhibits 
an excellent agreement with simulations in 2, 3, 4, and 
5 dimensions, nrovided that an accurate g{<7) is used as 
input .BBlla-ElEl On the other hand, if one is directly 
interested in obtaining reliable contact values gij(o'y ), 
then Eq. ( pl| ) is too crude. The simplest functional form 
of G that complies with ([^)-([To|) is a quadratic function 
of z: 



Girj, z) = Go(77) + Gi(77)2 + G2{ii)z\ 
where the coefficients are explicitly given by 

1 



Go(^) = 



1-77' 



Gi(7;) = (2-2^-277)5(^) 



2-r;/2 

1-7, ' 



1-77/2 



1 



(1 



■)ci-2 



77 



(12) 



(13a) 



(13b) 



(13c) 



In the one-dimensional case, Eqs. (13b) and (13c) lead 
to Gi = G2 = and we recover the exact result. For 
three-dimensional systems, if the SPT value is used for 
the one-component contact value, (7sPT(f) = (1 — 77/2 -|- 
77^/4)7(1 — 7,)^, we reobtain the SPT expressions for 
the mixture, cf. Eq. (^. On the othci; hand, if the 
much more accurate Carnahan-Starling^il (CS) expres- 
sion gcs{<y) = (1 ^ '7/2)/(l — is used as input, we 
arrive at the following expression: 



G(7/,2) = 



1 



1 3 77(1-77/3) 77^(1 -^/2) 2 

7; 2 (l-7/)2 ^ (l-7/)3 



(d = 3), 
(14) 

which is different from the BGHLL one and improves the 
latter for 2 > 1, as comparison with computer simula- 
tions will show. It should be noted that if one consid- 
ers a binary mixture in the infinite solute dilution limit, 
namely xi 0, so that 212 — > 2/(1 -I- (72/(^1), Eq. (1^) 
yields the same result foij-ai2(cri2) as the one proposed by 
Matyushov and Ladanyillj for this quantity on the basis 
of exact geometrical relations. However, the extension 
that the same authors propose when there is a nonva- 
nishing solute concentration, i.e. for xi ^ Q [cf. Eq. 
(19) in Ref. |l|l, is different from Eq. (|l|). We wiU come 
back to this point later when we assess the merits of both 
proposals. 

Of course, the quadratic form (0)is not the only 
choice compatible with conditions (|q)-(|lO|) . Another sim- 
ple possibility is to assume a rational function of the form 



G(77,2) 



l+[^l(7;)-l]^ 



(15) 
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Imposing Eqs. (|^)-(|To|), we get 

Boiv) - 1 - ^, 



2 l-(l-2''-ir;).g(CT)' 



(16a) 



(16b) 



1 2 - 7? - (1 ^ 77) (2 - 2''-^) .g(a) 
''^^'^^^-2 1- (1-2^-1,) g(a) ■ 

Other functional forms for G{r], z) coniplying with 
Eqs. (|)-(|lOl) are possible, but the choices (^ and (|l|) 
are sufficiently representative, so we will restrict ourselves 
to them in this paper. 




III. COMPARISON WITH SIMULATION DATA 

A. Contact values of the radial distribution 
functions 

In order to assess the utility of the new proposals for 
the contact values of the rdf, in Figs. ^]-^ we present re- 
sults for hard-core mixtures in d = 2, 3, 4, and 5, and 
the available computer simulation data. Figure ^ shows 
gijipij) as a function of for d — 2 and 77 = 0.6, accord- 
ing to the linear approximation (|l^), the quadratic ap- 
proximation (12), and the rational approximation (15). 



In the three cases we have used for g{a) the value ob- 
tained toni the Levin [6] approximant of Erpenbeck and 
Luban.Ej The only tabulated simulation data for gij{(Tij) 
in the case of hard disks that we are aware of are those of 
Ref. |3^. Hence, we have included in Fig. |l]the simulation 
results for the most asymmetric mixtures considered in 
Ref. |3^, namely 02! (J\ = ^ with x\ = 0.25, 0.5, and 0.75, 
and also simulation data extracted from Fig. 2 of Ref. 
|3^ . We observe that the quadratic and rational approxi- 
mations, both consistent with condition (10), are hardly 
distinguishable. The three theoretical curves practically 
coincide in the range of values of spanned by the sim- 
ulations. It would be interesting to carry out simulations 
extending to the region > 2 to assess the reliability of 
Eqs. (|l|) and (|l|). 

A comparison between theoretical predictions and sim- 
ulation values for three-dimensional mixtures is shown in 
Figs 



To carry out the computations in Eqs. (11), 
(|12|), and ( |15|) , we have used the CS contact value gcs(f). 
Figures |^ and ^ show that the universality assumption 
(^ is well supported by simulation data. Since the de- 
pendence of gij {oij ) on Zy is nonlinear (note that the 
curvature is different from that of the two-dimensional 
case), Eq. ([l^) only captures some kind of 'average' be- 
havior. Among the three quadratic functions, namely the 
SPT, the BGHLL, and Eq. Jl|), the best global agree- 
ment is presented by Eq. (|l^)7 The SPT prescription, 
Eq. (|^) with A = |, is consistent with condition (10), 



FIG. 1: Plot of the contact value gijipij) as a function of the 
parameter Zij — [uiUj / (Jij){a) / (a^) for hard disks (d = 2) at a 
packing fraction r) = 0.6. Tke open circles are simulation data 
for three binary mixturestJ with 02/ ci = \ and xi = 0.25, 
0.5, and 0.75. 'fke filled circles are simulation data for two 
binary mixture£j with aijox = | and x\ = 0.483 and with 
o"2/o"i = i and x\ = 0.219. The lines correspond to Eq. ( |ll| ) 
(dashed line), Eq. (^ (solid line), and Eq. (|l|) (dotted line). 



but its performance is not very good because it is pinned 
at a too high value at 2 = 1 (one-component case). The 
BGHLL prescription, Eq. (0) with A = i, is excellent at 
z = 1 (CS value), does a very good job for < z < 1, but 
clearly underestimates the simulation data for z > 1, as 
expected from the fact that the BGHLL is inconsistent 
with ( [l0| ) at z = 2. Our recipe ( |l^ ) is only slightly worse 
than the BGHLL for z < 1 but improves it significantly 
for z > 1. Finally, the rational function ( p^ is practically 
indistinguishable from the BGHLL for z < 1, is reason- 
ably good for 1 < z < 2, and is the best one in the case 
of the large-large rdf for disparate mixtures, as shown by 
Fig. H in the region z « 4. Of course, none of these ap- 
proximations is expected to be good enough ip-.thc .limit 
of extremely large values of z, where InG ~ z.qEjH The 
latter behavior could be incorporated by choosing an ad- 
equate functional form for G(?7, z) consistent with condi- 
tions (|^)-(|To|), but this does not seem to be necessary in 
the range < z < 4. Since the best global agreement in 
the range < z < 2 is provided by the polynomial func- 
tion (^^ , which has a structure similar to the well-known 
BGHLL prescription in the case of d = 3 [cf. Eq. (p^j, we 
favor its use, except perhaps for very disparate mixtures, 
where the rational approximation ([l^) is preferable. 

The refinements of the BGHLL|^exru:£ssions recently 
proposed lw-|Henderson and ChanQaoO and by Barrio 
and Solanao are not shown in Figs. ^ and |^ because they 
do not belong to the class of approximations satisfying 
the universality assumption (^. In addition, they are re- 
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FIG. 2: Plot of the contact value gij{o'ij) as a function of the 
parameter Zij = {aiaj /aij){a^) / {a"^) for hard spheres {d = 3) 
at a packing fraction r; = 0.49. Trhie symbols are simula- 
tion datarfor the single fluid (circle) ,t3 three binary mixtures 
(squares)Ea with 0-2/0-1 = 0.3 and xi =p-j0.0625, 0.125, and 
0.25, and a ternary mixture (triangles )Ej with 0-2/0-1 — |, 
cs/cri = I and xi — 0.1, X2 ~ 0.2. The lines are, from bot- 
tom to top at the right end, Eq. (0) ( ), BGHLL (-■-•). 

Eq. ^ (• ■ •), Eq. (|l|) (-), and SPT ( ). 



stricted to the case of binary mixtures. Both approxima- 
tions differ in practice from the BGHLL only in gii(a-i), 
-where species 1 refer to the big spheres (cti > 172). Fig- 
ure ^ sho-ws 511 (ci) versus ij for the three-dimensional 
binary mixture xi — 0.005, 02! (J\ = ^ (-which corre- 
sponds to zii = 3.457). The figure confirms that the best 
agreement is obtained -with the rational approximation 
(p^). Henderson and Chan's approximation, -which incor- 
porates the exact asymptotic behavior Ingii(cri) cti, 
gives too high values. Barrio and Solana's expression im- 
proves the BGHLL value, but is slightly -worse than the 
quadratic approximation (|l2|). 

Figure |5| presents a plot of (j'i2(cri2) as a function of 
r] for the three-dimensional binary mixture characterized 
by Xi = 0.00415 and o^jox = \ (-which corresponds to 
a value of Zi2 — 1.210) as given by different approxi- 
mations. The rational approximation given by Eq. (|l5| ) 
has not been included in the figure to avoid overcro-wding 
of the curves, but it is practically indistinguishable from 
the BGHLL approximation in this case. Clearly the best 
agreement bet-ween theory and simulatioB-is provided by 
the approxixnations of Barrio and Solana,cll of Henderson 
and Chan,Ej and by our Eq. (p^), -which are all of com- 
parable accuracy and certainly superior to the approxi- 
mation proposed in Eq. (19) of the paper by Matyushov 
and LadanyiO 

Of course, the most physically relevant situations cor- 
respond to hard spheres [d = 3) and, to a lesser extent, 
disks (d = 2). On the other hand, it seems desirable 



FIG. 3: Plot of the contact value gij{(Tij) as a function of the 
parameter Zij = {ataj /aij){a^)/{a^) for hard spheres {d = 3) 
at a packing fraction rj = 0.314. Trhe symbols are simula- 
tion data for three binary mixturedij with 02/01 = i and 
xi = 0.00311, 0.00415, and 0.005. The lines are, from bottom 

to top at the right end, Eq. (^ (- - -), BGHLL (- ), 

Eq. (|l|) (-), SPT ( ), and Eq. (|l|) (■ ■ •)• 

that a proposal for gij [uij ) be valid for any dimension- 
ality d. Moreover, a number pf. repent.. papers -deal -with 
systems of hard hyperspheresrjO'LaSHcilo'H Figures 
sho-w 512(0-12) versus rj for binary mixtures in d = 4 
and d — b. The contact values of the one-component sys- 
tem that -we have used in the computations have been ob- 
tained from the EOS derived for these systems by Luban 
and Michels.E3 The values of the parameter Z12 are in 
the range 0.4 < Z12 < 1 for the cases considered in Figs. 

It is observed that in this range the quadratic ap- 
proximation dl^ ) exhibits an excellent agreement -with 
the simulation data. 



B. Equation of state 

Having examined the accuracy of the proposed contact 
values, -we -will no-w consider their performance in terms 
of the compressibility factor they lead to. In this regard, 
Eq. (jlj) has the advantage over Eq. ([l^) that, -when in- 
serted into Eq. (|l|), one gets a closed expression for the 
compressibility factor, in terms of the packing fraction 77 
and the first fe-w moments (cr"), n < d. This expression 
is meaningful even for polydisperse mixtures. The result 
is 

Z„,(ry) - 1 + 2^^-^^ [2(5o - 25i + S'2) + (5i - 52)ry] 
1 -?7 

+ [Zsiv) - 1] [25i -S2 + 2''-^{S2 - 51)77] ,(17) 

-where Zs{rf) — l-\-2'^~^rjg{a) is the compressibility factor 
of the one-component system and the coefficients Sm are 
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FIG. 4: Plot of the contact value gii{ai) as a function of the 
packing fraction 77 for the three-dimensional binary mixture 
xi — 0.005, pOta/o-i — i (211 = 3.457). The symbols are simu- 
lation data.LJ The lines are, from bottom to top at right 

end. Eg. (0) ( ), BGHLL {- ■ - j), Barrio-SolanaH (■■■), 

Eq. (Qi^jri, SPT ( ), Eq. @ (■ ■ •), and Henderson- 
Chan 




FIG. 5: Plot of the contact value t;i2(o-i2) as a function of the 
packing fraction 77 for the three-dimensional binary mixture 
XI = 0.00415, (Wo-i = i {zi2 = 1.210). The symbols are 
simulation data.cJ The lines are, from bottom to top at the 

right end, Eq^Jl^) ( ), BGHLL (- rrnr^jrfq- & (-)- 

Barrio-Solanan^--), Henderson-rehanHElBM 1 
(- •■ - •■), and Matyushov-LadanyitJ (■ ■ ■). 




SPT 



given by 



-(d-m) (^'^ ^ ™] (ct"+"')(ct'^^") 



In the two-dimensional case, Eq. ( |l7| ) becomes 

1 



(18) 



Zr-niv) = T— + 7-1 

1 — 77 (ct"^ 



,2 r 



1-77 



(d = 2). (19) 



It is worth noticing that this EOS coincides with the 
one obtained from Eq. (|Tl|) for d = 2. This illustrates 
the fact that two different proposals for the contact val- 
ues Qiiicij) can yield the same EOS when inserted into 
Eq. Let us analyze this point with more detail. Sub- 
tracting Eqs. (0) and ([l^), one has 



AG(77,z) 



1-77/2 



1-77 



(1 



z(l-z), (20) 



where AG'(77, z) denotes the difference between the linear 
and the quadratic approximations. Thus, the compress- 
ibility factors obtained from the linear and quadratic 
forms for G(77, z) only differ by a term proportional to 
J2i j ^i^j'^ij^iji^ ~ ^ij)- It turns out that this term 
vanishes in the two-dimensional case, so the linear and 
quadratic approximations yield the same EOS (|l9|). This 
fact also shows that a rather crude approximatiQnj.suab. 
as Eq. (pTl) may lead to an extremely good E0S.Noe3'El 




FIG. 6: Plot of the contact value (?i2(o"i2) as a function of 
the packing fraction r; for the four-dimensional binary mix- 
ture x\ — 0.5, (fS|/cri = i (212 = 0.706). The symbols are 
simulation data.E3 The lines correspond to Eq. (0) (dashed 
line), Eq. (|l|) (dotted line), and Eq. (|l|) (solid line). 



For three-dimensional mixtures, Eq. ( [17] ) becomes 



^m(?7) = 



1-77 



1-77 



1 



\ — rj 



(d = 3). 



(21) 
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FIG. 7: Plot of the contact value (7i2(cri2) as a function of 
the packing fraction rj for the four-dimensional binary mix- 
ture a;i = |, cr2|/«ri = 0.25 {zi2 = 0.402). The symbols are 
simulation data.tj The lines correspond to Eq. (^l|) (dashed 
line), Eq. (|l|) (dotted line), and Eq. ^ (solid line). 




FIG. 8: Plot of the contact value gi2{o'i2) as a function of 
the packing fraction rj for the five-dimensional binary mix- 
ture a;i = i, (J2pZfiri = 0.75 (zi2 = 0.912). The symbols are 
simulation data.a The lines correspond to Eq. (^l|) (dashed 
line), Eq. (|l|) (dotted line), and Eq. (Ili) (solid line). 



In particular, when the CS EOS Zs{r]) = {1 + i] + t]'^ - 
r/^)/(l — ry)'^ is used as input, we get the following ex- 
tended CS EOS: 



(1 - 77)2((73)^ 



(22) 

where the compressibility factor associated with the BM- 



FIG. 9: Plot of the contact value t;i2(o"i2) as a function of 
the packing fraction 77 for the five-dimensional binary mix- 
ture a;i = i, (T*/<Ti — 0.5 {zi2 = 0.687). The symbols are 
simulation data.u The lines correspond to Eq. (^l|) (dashed 
line), Eq. (|l|) (dotted line), and Eq. (|l|) (solid line). 



CSL EOsli'i is given in the present notation by 



Z 



BMCSL 



l-T] (l-77)2(a3) (l-,y)3(cr3)2 " 

(23) 

Note that Eq. (|2|) is different from the extended CS EOS 
obtained from Eq. (11), namelyMl 



■^ccs-i(??) = ^bmcsl(?7)^ 



773(^2) 



(l-77)3(a3)- 



(24) 

Since simulation data indicate that the BMCSL EOS 
tends to underestimate the compressibility factor, it is 
obvious that the performance of Zccs-i is, paradoxically, 
better than that of .^cCS-ii- This is explicitly shown in 
Fig. |l^ where, for comparison, we have also included the 
theoretical results that follow from tho, recent (very ac- 
curate) proposal by Barrio and SolanaEj for the EOS. So 
once again we find that better contact values do not nec- 
essarily lead to better compressibility factors, as already 
seen in the two-dimensional case. One plausible explana- 
tion for the better performance of .^eCS-i with respect to 
■^oCS-ii might reside on the use of the CS EOS Zs(77\yJDut 
we have checked that if the Pade [4,3] of SanchezCJ or 
the very accurate EOS of Malijevsky and VeverkaEa are 
used instead, the results do not change. Thus the source 
of this effect is a fortunate compensation of errors in the 
linear approximation ( pi] ) related to the fact that the 
compressibility factor involves a weighted average of the 
individual contact values gij {aij ) [cf . Eq. (^] . This argu- 
ment is suggested by the following observation. In Fig. |ll| 
we present a plot of the ratio G{rj, 2:)/Gbghll(?7, z) as a 
function of z for 7, = 0.49 with G(r], z) given by Eqs. 
( pT| ) and (|l^). For completeness, also the ratios corre- 
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spending to the SPT (|^) and the rational approxima- 
tion ( |l^ ) are plotted. Properly reduced simulation results 
have also been included. From this plot it is fair to con- 
clude that the approximation given by Eq. ( |l2|) is globally 
more accurate than those obtained with cither Eq. 
or with Gbghll(?7, z) and that the BGHLL contact val- 
ues are better than the linear approximation. However, 
the quadratic approximation always underestimates the 
simulation results (and therefore it is only natural that it 
will always produce an underestimation of the compress- 
ibility factor), while both the linear approximation and 
the BGHLL approximation overestimate the simulation 
results if < 1 but underestimate them if > 1. The 
net result is that ZeCS-ii(??) is in poorer agreement with 
the simulation results for the compressibility factor than 
either ^oCS-i('7) or 2'bmcsl('?), the extended CS EOS 
obtained from Eq. (^l|) providing the best overall agree- 
ment. Nevertheless, an important asset of .^cCS-iii not 
shared by either Z^cs-i or ^bmcslj is that it predicts 
demixing. This result provides further support to the 
analysis performed by Regnaut, Dyan, and Amokranecj 
in which the verification of condition is of key impor- 
tance for the existence of demixing. We will address this 
issue in more detail elsewhere. 




FIG. 10: Deviation of the compressibility factor from the 
BMCSL value, as a function of the packing fraction ry for 
an equimolar t^ree-dimenpiffinal binary mixture with a^jax = 
0.6. The openEl and fillecO circles are simulation data. The 
lines are, from bottom to top at the right end, the eCS-II 
EOS ( p2[ ) (- — ), the EOS obtained from the rational approx- 
imation (|l|) (■ ■ ■), ther«CS-I EOS ( ), and Barrio 

and Solana's EOS (--OM 

For d = 4 and d = 5 (not shown), however, the com- 
pressibility factoraJcrivcd from the linear approximation 
given by Eq. ([ll|)ll3 turn out to be slightly less accurate 
than those obtained from the use of either Eq. (|l|) or 
Eq. (15). Thus, for these high dimensionalities, the ex- 
pectation that a better overall behavior of the contact 
values correlates with a better performance of the asso- 
ciated EOS actually holds true. 




FIG. 11: Plot of the ratio G(?7, z)/Gbghll(»7, z) as a func- 
tion of the parameter z for hard spheres id = 3) at a pack- 
ing fraction r\ — 0.49. Hie symbols are simulation data fefl 
the single fluid (circle) ,c3 three binary mixtures (squares)E3 
with (T2/o"i ~ 0.3 and xi 0.0625, 0.125, and 0.25, and a 
ternary mixture (triangles)E3 with (72/cri — |, crx /cn = ^ and 

xi — 0.1, X2 = 0.2. The lines correspond to Eq. ( |ll[ ) ( ), 

Eq. (|l|) (-), Eq. (0 (■ • ■), and SPT ( ). 



IV. CONCLUSION 

In summary, in this paper we have introduced a univer- 
sality assumption, Eq. (|^), for the contact values of the 
rdf of a hard-sphere mixture with arbitrary number of 
components and arbitrary dimensionality. Three known 
consistency conditions, Eqs. (||)-(|4|), allow us to fix the 
values of the universal function Gij], z) at z = [cf. Eq. 
(|)], z = 1 [cf. Eq. (|)], and z = 2 [cf. Eq. (0], the two 
latter in terms of the contact value of the one-component 
rdf. This implies that any reasonable three-parameter 
form of G{r], z) as a function of z can provide a very 
good approximate representation of gij (uij ) regardless of 
the number of components, provided a good EOS for the 
single fluid system is used. We have illustrated this pos- 
sibility with two specific proposals: a quadratic function, 
Eq. (|l^), and a rational function, Eq. (|l5|). In d = 1, 
they reduce to the exact result, while for d = 3 they rep- 
resent an improvement over the BGHLL values, as well as 
over those of their refinements, in the quantitative agree- 
ment with the simulation results. For d = 2, 4, and 5, 
they compare rather well with the (few) available simu- 
lation results. Their potential use in connection with the 
generation of the entire rdf gij{r) for three-dimensional 
mixtures within the Rational Function Approximation 
method is currently under investigation. 

The relationship between the thermodynamic proper- 
ties and the contact values of the rdf in hard-core fluids 
is in principle so straightforward that the importance of 
explicit and accurate approximations for the latter can 
be hardly overemphasized. In one-component hard-core 
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systems it is certainly true that a more accurate contact 
value of the rdf leads directly to a better EOS. On the 
other hand, our results in this paper indicate that, due to 
the fact that for mixtures the EOS (or equivalently the 
compressibility factor) involves a summation over species 
indices in which the different contact values are included, 
one does not always obtain a more accurate EOS from 
seemingly better approximations to gij{(Jij). In fact, as 
exemplified in the case of hard-disk mixtures [d = 2) , it 
is possible to obtain exactly the same EOS with two dif- 
ferent approximations for the contact values of the rdf. 
Further, the poorer agreement of .Zecs-ii with simulation 
data than either Z^cs-i or 2'bmcsl mentioned in the pre- 
vious Section is also a reflection of the above assertion, 
the reason resting on a 'fortunate' compensation of er- 
rors. In any event, in some specific applications (e.g. the 
Enskog kinetic theory) it is only the contact values of the 
rdf that are required. In this respect, it is fair to con- 
clude that our two new proposals provide in general a 
reasonably accurate approximation to gij{(Tij) (as com- 



pared to the available simulation data) for a hard-core 
mixture with an arbitrary number of components and 
arbitrary dimensionality. Of course the scarcity of simu- 
lation results for these systems precludes a more definite 
conclusion. In any case, we hope that the availability 
of the new (explicit) expressions for the contact values 
of the rdf of hard-core mixtures in any dimensionality 
may serve as a further motivation to carry out yet more 
simulations of these systems. 
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